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ABSTRACT 


Using present linear theory results for a rigid simple 
plenum fluid-suspension system, it was possible to determine 
the accuracy of a non-linear model and of analog computer 
Simulation when using small perturbation amplitudes. 


By increasing tne perturbation amplitudes it was possible 
mo Snow tnat the range of validity of the linear theory solution 
increased as the lead time constant increased and also as the 
lead-to-lag time constant ratio increased. 


furtner increasing of the perturbation amplitude snowed that 
tne maximun allowable perturbation amplitude increased as the lead 
time constant increased and also as the lead-to-lag time constant 
ratio increased. 


Increasing the perturbation amplitude beyond the range of 
validity of tne linear theory snowed that tne peak dimensionless 
acceleration and the peak dimersionless change in vehicle- 
Suspension~guideway clearance decres#sed. 


Tnesis Suvervisor: H. H. Richardson 
Associate Professor of Mechanical bngineering 
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Denotes time derivative 


— Denotes dimensionless value obtained by dividing by the 
equilibrium valve (unless noted above) 
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INTAODUCTION 


emer eee 


Phrovusnout his history, man has continually oti aca ome 
himself and his soods from one place to anotner Witn ever- 
increasing speed. The horse, Snip, railroad train, automobile, 
ano airvlane all offered increases in speed and were, therefore, 
readily accevted. 

We have now, however, reached a paradox in the field of 
short haul (up to distances of four to six hundred miles) 
transportation. The fastest modetof transportation 1S often in 
reality the most time consuming. This is the result of congestion 
at airports, and the advent of the "“super-jets" and "“air-busses" 
Will only compound the problem. An alternate solution involving 
smaller, but more frequent, aircraft merely moves the congestion 
from the terminal building to the runway, also increasing costs 
Sewa result of the decreased payload to gross welght ratio. 

Efforts to increase the speed of the surface transportation 
modes have thus been attempted in many areas or tne world. Inese 
efforts have, however, pointed up two major limitations: 
mechanical suspension deficiencies and hign power requirements, 

The suspension deficiencies are of prime importance as these, 
Coupled with road-bed irregularities, produce intolerable 
accelerations and heave motions at speeds of two to three hundred 


miles per hour. Efforts to reduce the road~vded irregularities have 








proven to be both expensive and of short duration because of 
uneven loading and settling. Also, the limit appears to have 
been reached in the design of mechanical suspension systems 
with no major break-tnroughs expected in the near future. 

With these limitations in mind, it was decided to attempt 
a fresh approach to the suspension system for a hich-speed 
Sround transportation vehicle. The result of this new approach 
Was the fluid-suspended vehicle. A form of this type of 
suspension is shown schematically in figure one (1). As shown, 
the vehicle does not make physical contact with the Suidevay, 
bul ratner, floats on a cushion of compressed air. It is this 
cushion of air whicn gives the suspension its dynamic 
Mearacteristics. This cushion of air acts as a spring with its 
corresponding lead time constant and also acts to impart a 
lags time constant, as will be snown later. 

This fluid-suspension system suggests advantages over the 
mecnanical suspension in the consequent financial savings in 
meigaeway construction. Through the use of low pressure air, 
of approximately one pound per square inch, gause, the high unit 
loadings encountered in tne mecnanical suspension can be 
eliminated, thus prolonging guideway life. Also, locally 
rousn materials, such as concrete, way be used as the guidewa 
material since there are no wheels or tires to Wear out. 


One disadvantage of the mechanical suspension cited 
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above is the hign power requirement necessary to overcome 

wind drag on the vehicle. This power requirement is still 
present witn the use of a fluid suspension and. 1S ena i. 
augmented by the power required to support the vehicle. At 

the hign overating speeds envisioned, however, this hovering 
power is expected to be small compared With the’total power 
requirements. At these speeds the weight of the vwelbicle coud 
conceivably be supported by compressed air provided either 

by aerodynamic means or by induction of air from the stagnation 
region at the front of the vehicle. In this way the hover-power 
requirements could be reduced and thus the total potvrer requirement 
would approach that of a mechanically-suspended vehicle. 

Because of the anticipated advanteges of a fluid-suspended 
eiicle, a linearized model was formulated to further explore 
the vehicular dynamics. This exploration was begun by striving 
tO combine tne basic equations describing tne suspension. 
these equations included ones describing flow Ante And "ong o. 
moiewcwsnion, the continuity eguation for the cushion volume, 
Newton's Second Law, and the equation of state of tne fluid 
and are the same as those used for the non-linear model in the 
Detailed Procedure. The simplest transfer function which can 
be formed from these equations, assuming the venicle-suspension~ 
gsuideway clearance to be equal around tne periphery of the 


fmeiton, is 











a (1) 


where, for tne simple risid plenim suspension, 


NS ee (2) 


is the incremental cnange in suspension force, 


Ah = h-h 3) 


is the incremental change in the vehicle~suspension~guideway 


clearance, 


a 8 é 


: e 1+1/6F 


(4) 
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is the static stiffness of the suspension, 


t= (5) 


is the lead time constant and physically represents the time a 


fluid particle spends in the active cushion volume, 





(6) 


is the lag time constent and pnysically represents the charging 
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las of tne entire cushion volume considered as a fluid 


Capacitance, 


oe eee 
Se (7) 
Poe Pa 


is the ratio of maximum incrementel force Which the suspension 


Gan carry to the equilibrium force, 
so aes 
ee (8) 


is the dimensionless meximum active volume at equilibrium, and 


ae W 
W= <a ( 9 ) 
. tae ace € 





is the reciprocal of the time a fluid particle spends in the 
entire cushion volume at equilibrium. 

i[aeom dave. on the distribution of guideway irregularities, 
both in magnitude and average distance between these perturbations, 
it is possible to describe the perturbations as a sinusoidal 
Gisturbance. This disturbance would nave a frequency equal 
to the venicle's forward velocity divided by tne averace distance 
between perturbations. This wovld suggest that the transfer 
function should be solved for sinusoidal disturtances. when 


this is done, it yields 
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FIGURE 2, Fluid-Suspension Vehicle Moving at Velocity "v" 
over a Guideway with Averase Distance between 
Perumroava one '):! 
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where 

/A = the dimensionless vehicle acceleration 
lH = the dimensionless change in vehicle-suspension-~ 


suldeway clearance 


aa = ta (12) 
is the dimensionless disturbding frequency, 
mre Or (13) 
Mestre dimensionless ee time constant, 
re == Soe (44) 


1S the dimensionless lag time constant, and 


= [ae pe 
a, = [es = om | bee Gig) 
a eee, r/o 








is the natural frequency of the suspension. 


The solution of the dimensionless acceleration shows that 
peak values occur at or near the natural frequency of the venicle 
and at very large frequencies when the dimensionless lead time 
constant is less than one. By equating these two peek values, 


an optimum lag time constant can be obtained as a function of 


the lead time constant. 





fe 
1+ Ve a 
Ine at = Ht (16) 
Pp Ao a o> oo 
2 
ye 
yi = 1 (17) 
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Using, this optimum value of the lag time constant, it is 
How possible to solve equations ten (10) and eleven (14) for 
the peax values of the dimensionless vehicle acceleration and 
dimensionless change in vehicle~suspension-suidevay clearance 
fOr a Darticular lead time constant. These solutions are 
shown graphically in figure three (3). Using this figure and 
data on the distribution of guidevay irresularities, both in 


magnitude and in average distance between tnese perturbations, 


it is now possible to design a fluid-suspension system and 
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mebrermine its dynamic characteristics. 


ms. a result of Limitations of linear Solutience ic ever 
questions arose as to the range of validity of the solutions 
and also as to the response of the vehicle beyond this range. 
It is the purpose of this thesis to answer these ‘questions for 
the case of a simple rigid plenum fluid-suspended vehicle and. 
to investigate the condition where the lag time constant equals 
Bhe lead time constant. It is at this point that the linearized 


transfer function becomes infinite. 
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By using a simple rigid plenum fluid-suspension systen, a 
non-linear mathematical model was formulated and simulated on an 
analog computer. The details of this formulation and simulation 
are contained ‘Via ‘iced Lis) PROCEDURE. Daas Simulation was 
then used to determine the maximum acceleration and the maximum 
change in vehicle-suspension-guideway clearance for a range of 
lead time constants and their associated optimum lag time 
constants. By checking these values for small perturbations, 
the simulation was verified. The perturbation amplitude was 
then increased until the maximum value allowed by the system 
was aeined. This occurred either when the vehicle struck the 
guideway or When the cushion pressure equaled or exceeded the 
source pressure. 

Ine simulation was then used to determine the maximum 
acceleration and the maximum change in vehicle-suspension- 
Suideway clearance for the condition of a fixed lead time constant 
with values of the lag time constant avproaching thet of the 


mesq time constant. 








RESULTS AND DISCUSSION 
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rigures four (4) and five (5), which contain the results of 
varying the lead time constant with the optimum lag time 
eonstant, show that the peak dimensionless acceleration and the 
peak dimensionless change in vehicle~suspension-guideway 
Clearance obtained for a dimensionless perturbation height 
of 0.0] correspond very closely to the linear theory solution 
for all values of the lead time constant. These variances are 
probably the result of round-off errors in setting up the problem 
on the computer. 

These figures also show tnat as the lead time constant is 
increased, the range of validity of the linear solution increases 
to a dimensionless perturbation height of 0.2, and tne maximum 
allowable dimensionless perturbation height increases to C.4. 

it is also noted tnat the peak dimensionless acceleration 
and the peak dimensionless change in vehicle-suspension-guideway 
clearance decrease as the dimensionless perturbation heignt is 
increased above tne linear solution range. 

Fisures six (6) and seven (7), which contain the results 
of allowing the lag time Jemeiv ane to approach and equal the 
lead time constant, snow that the peak dimensionless acceleration 
and peak dimensionless change in vehicle~-suspension-sguidevay 


Clearance obtained for a dimensionless perturbation height of 








Ox 


See Sp 


Y5* 0.01 
0-02 


O.07 KR 


PEAK DIMENSIONLESS ACCELERATION 


~~? 


jo iO 





DIMENSIONLESS LEAD TINE CONSTANT 
wt ® 

(bg 
4 


FTGURE 4, Peak Dimensionless Vertical Acceleration Versus Lead 
Time Constant and Dimensionless Perturcation Height 
for Optimum Lag Time Constant 
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PEAK DIMENSIONLESS CHANGE IN VEHIC 


0.4 — 


(een ees ee 
om 107 / 


-. 


Dita NSiOtieess LEAD Gils CONSTANT 


FIGURE 5, Peak Dimensionless Change in Venicle-Suspension- 
Guidewey Clearance Versus Lead Time Constant and 
Dimensionless Perturbation Heisht for Optimum Lag 
Time Constent 
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FIGURE 6, Peak Dimensionless Vertical Acceleration Versus 
Lead-to-Lag Time Constant Ratio and Dimensionless 
Perturbation Height 
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FIGURE 7, Peak Dimensionless Change in Vehicle-Suspdension-~ 
Th 4 
Guidevay Clearance Versus Lead-to-Laz Time Constant 
Ratio and Dimensionless Perturbation Height 
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mol correspond very closely to the linear solution. There is 
one exception, however, in the case where the lead and las time 
constants are equal. At this point linear theory predicts 
infinite values. However, the computer showed a finite value 
for a dimensionless perturbation height of 0.01. 

When a disturbance was introduced and then removed, the 
system continued to oscillate indefinitely with a peak 
acceleration of 15.36 feet per second saquared and a peak change 
in vehicle-suspension-guidevay clearance of 44 per cent of the 
eavilibrium height. 

These observations sugsest that the linear theory solution 
is only valid for infinitely small disturbances. 

These figures also show that as the lead-to-lag time 
constant ratio increases to two, the range of validity of the 
linear solution increases to a dimensionless perturbation 
height of 0.2 and the maximum dimensionless perturbation heignt 
allowed increases to 0.2. 

As elso noted in the other case, the peak dimensionless 
acceleration and the peak dimensionless change in vehicle- 
Suspension-suideway clearance decreases as tne dimensionless 
perturbation neignt is increased above tne lineer solution ranse. 

ithe fact that tne peak dimensionless acceleration and the 
peak dimensionless change in vehicle-susvension-suidetiay 


clearance decreases for dimensionless perturbation heights 








beyond the range of validity of the linear theory solution 


implies tnat there is a bias in the system. This blas can be 


shown quite clearly in the case where the lag time constant 


equaled the lead time constant and the system oscillated 


indefinitely at zero perturbation heisht. 
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CONCLUSIONS AND RECOMMENDATIONS 


fi. When the lead time constant is increased using the optimum 


lag time constant 


a. The range of validity of the linear solution increases 
iow The maximum allowable dimensionless perturbation height 
increases 


a When the lag time constant is allowed to approach the lead 
time constant 
QB. The range of validity of the linear solution decreases 
uO Zero 
or The maximum allowable dimensionless perturbation height 
decreases 
B. When the dimensionless perturvation height is increased 
beyond the range of validity of the linear theory, the 
peak dimensionless acceleration and the peak dinensionless 


change in vehicle-suspension-guideway clearance decrease. 


it is recommended that, wnen the new Mechanical Engineering 
Devartment Analog Computer facilities become available, this 
model be used as a basis for investigation of the bias noted 


in the system. 
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DETAILED PROCEDURE 


DEVELOP ANT OF NON-DINSAR MODS. 


Samia 


Tne continuity equation for the cushion volume can be 


eeeacGecen as 


We = Povo Bae + Vo | (18) 


Assuming that tne density perturbations are small, this can 


be nondimensionlized about equilibrium as 


¢ @ 


~~ —— 8 oe 


& Ne ~ WwW) we -- Vo fe (us)) 


wnere 


7 Wi 
e 
Pee Ce 
Since a practical suspension would Keep the pressure 
differences as small as possible, the incompressiole flow 


equation can be used to describe the flows into and ovt of 


tne cushion volune as 


SP pow fe OO = ee ee ee 


Nig = Ere 21o OD ah v9) 
Wo= WLC 5/205 9 (PoP, ) (215) 





At equilibrium tnese flow equations become 





We = Se = A Cael2 Pio (Pe-Poe) (22) 
Se eee cay) ce Poep mau (23) 


The flow equations can be nondimensionalized by dividing by 


the equilibrium conditions. 


-- S s -c 
Wo=— = |= (24) 
: ws Ps Poe 
OR nd 
—- WwW c ta 
Woe on = a= 25) 
e ce “a 


Newton's Second Law can be used to equate the cushion 


pressure and tne acceleretion of the vehicle as 


F=A(v +-0_) = m( sty.) (26) 


Nondimensionalizings about eauilibrium oroduces 
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Poe Pa Te es 428 
where 
pe 2 : 
Y2 = he (29) 


Taking the time derivative of the above vields 


DB, AS es 
- =~: ¥ (30) 
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Combining the above witn tne flow equations to eliminate 


the pressure terms produces 








re = (34) 

a Hd 

VW f= ; es = Yo (32) 
wnere 

= ate 

ES aa (7) 
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The cushion volume can be written as 


V 
Cc 


where 


V+ Ah 
O 


Ce 


is the cushion dead volume. 


35) 


(34) 


Dividing by the equilibrium cusnion volume to nondimensionalize 


gives 


where 


| 


sf 


Again assuming the density perturbations are small, and 





=e of =V, + Ah 
© GE Se 
Anh 
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iM) 28s 
i; h 
S ce 


a 


(25) 
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assuming a polytropic process, the equation of state of the 


Panid is 


Po Yo _ + 
D = Yeo (38) 


ce - Poe 


Combining with equation thirty (30), the above becomes 
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Substitutins the above expressions into equation nineteen 


(19), the finel form of the non-linear mathematical model is 
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DeVeLOPnih ! Or nk ANATOG COMPUTE SIMULALLION 


The analog computer used for this simulation consisted 
Of Philorick Universal Linear Operators, model K5-U, and 
Philbrick Multiplier-Dividers, model K5-M. The outputs were 
displayed on a EB paeenier Displey System, model 5934. 

The non-linear mathematical model can be solved for the 
hishest order derivative of the vehicular displacement in 


merms of the lower order derivatives as 





Poe Pa net 7 
YPoe 8 ¢ 
Sa: Vee ee 
- So = by 
a aca ae Ah | (414 ) 


Using this equation, the operational bice™: diagram can 
tnen be Duilt up by vsing integrators to calculate the lower 
order derivatives, trus closing the feedback loovs and completing 
maemesOlurion. Tne onerational block diagram thus formed, 
Pocovumims tne scaling factor is one volt equals one dimensionless 


4 


wees Shown in fisure eisnt (8). 


In 


vn 
In determining the voltase scaling faucor to ve used for 
the actual simulation, tnere were two major considerations 


resulting fron limitations of the analog computer components, 
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FIGURE &, Operetional Block Diagram of a Simple Higid Plenum 
Suspension Assuming a Unity Voltase Scaling Factor 








The first consideration involves the fact that the comvuter 
Gimensionless vehicle acceleration (vehicle acceleration divided 
meethe acceleration due to gravity) is Sens toe ve rather 
small and therefore near the noise level of the computer. 
iaeretore, the scaling factor should be as high as possible. 

The second consideration determines the upper elem, 
as the computer dimensionless clearance (vehicle clearance 
Meieht divided by the equilibrium clearance height) could 
dovole in magnitude. Thus the scaling factor cannot be too 
large as tnis value migne then so beyond the maximum operetional 
mevel of the computer which is fifty volts. Thus a compromise 
voltage scale of twenty-five volts per dimensionless unit was 
used. 


The time scaling used was real time. 





DeVELOPMENT OF DATA POINTS 


The computer Simulation was used to obtain two sets of 
Gata throuscn the use of two different sets of initial 
Bonditions. 

The first was through the use of the optinum Lag sbame. ae 
determined from figure three (3), for an assumed lead time 
constant. The remainder of the design parameters of this set 


of data could then be calculated assuming the following 


constraints: 

fle The ratio of tne maximum incremental force 
to the equilibrium force (8F) is one. 

Le The equilibrium clearance om ) iis one inch. 

oF The minimum allovrable dimensionless active 
volume (A) is 0.1. 

Mh, Only discrete equilibrium cushion pressures 
(Doo) Gime toe voncs. and sto. 7 Dsioeane 
used. 


Dividing equation five (5) by equation six (6) to obtain 
eS eo ome se gs | (2) 


it is possible to determine the lowest equilibrium cushion 


pressure which will satisfy the constraint on the dimensionless 





active volume and to then solve for tne dimensionless active 
volume usins the next larsest discrete value of the equilibrium 
cushion volume. 

The dimensionless dead volume (V.) was then calculated 


from 


ae ee (43) 


which can be obtained from equation thirty-five (35) at 
equilibriun. 


The natural frequency (w-,) can then de calculated fron 


a = 


ie) 








n 
Wo ( ity ) 
e eee 
Z4 
which can be obtained from combining equations five (5) and 
thirteen (13). 
Using these values, the computer was programed and 
tested against the linear tneory solution using a dimensionless 


perturbation height of 0.01. If the results corresponded with 


the linear tneory solution, tne dimensionless perturbation 





height was increased until eitner the venicle struck the 
suideway or the cushion pressure excected tne source pressure. 
The second set of data was Ape e in the same manner 
except that the lead to lag time constant ratio was assumed, 
thereby determining the lag time constant, rather than using 


the optimum les time constant as above. ; 





SUMMARY Of DATA 


the followings two tables contain the summary of tne data 


sd as the basis for ficures four (4) through seven (7). 
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TABLE 1, Summary of Data for Various Lead Time Constants 
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SAIPLS CALCULATION 
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wt - 0.4 (assumed ) 

> - 0.3 (from fisure two) 
o4 = O.4 ae Li, 0 
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Determinetion of the minimum equilibrium cushion pressure 
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Miererore, the equilibrium cushion pressure will be 15.7 psie, and 





“a a J) ae . 
ir ieee 8 OT oo es 2a 
OF fel ew Pa ee Oe 1 2 
Zh YPog 4a,/ep  O-2 (1-4)45.7) 144 
Vo = 4 -A=14 - 0.182 = 0.8418 





lege ae 
J At /bF 


ae AD 

ze 8 

Ww no AO a2) Nesom) = 8.94 sec at! 
e ia 0.4 


Spf nD ee A 
=f T1742) ./ T47.= 19-05 see 


From figure two (2), the peak dimensionless acceleration and 


the peak dimensionless change in vehicle-suspension-fuideway 


clearance are: 
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For a dimensionless disturbance of Yofh 


the ooserved values are: 
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